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ABSTRACT 

Dynamo action owing to helically forced turbulence and large-scale shear is studied using direct numerical 
simulations. The resulting magnetic field displays propagating wave-like behavior. This behavior can be 
modelled in terms of an aQ. dynamo. In most cases super-equipartition fields are generated. By varying the 
fraction of helicity of the turbulence the regeneration of poloidal fields via the helicity effect (corresponding to 
the or-effect) is regulated. The saturation level of the magnetic field in the numerical models is consistent with 
a linear dependence on the ratio of the fractional helicities of the small and large-scale fields, as predicted by 
a simple nonlinear mean-field model. As the magnetic Reynolds number (Re^) based on the wavenumber of 
the energy-carrying eddies is increased from 1 to 180, the cycle frequency of the large-scale field is found to 
decrease by a factor of about 6 in cases where the turbulence is fully helical. This is interpreted in terms of the 
turbulent magnetic diffusivity, which is found to be only weakly dependent on Re^- 

Subject headings: MHD - turbulence 



1. INTRODUCTION 

Several classes of turbulent astrophysical bodies including 
stars with outer convection zones and spiral galaxies can con- 
tain pronounced large-scale magnetic fields. The scales of 
these fields are much larger than the scale of the energy- 
containing eddies of the turbulence responsible for produc- 
ing these fields. Explaining such fields is an important aim 
of dynamo theory. Today, turbulence simulations begin to re- 
produce the phenomenon of large-scale field generation (e.g. 
Brandenburg 2001; Brandenburg et al. 2001, 2008a; Yousef 
et al. 2008a,b; Kapyla et al. 2008, 2009a; Hughes & Proc- 
tor 2009) and allow making contact with mean-field dynamo 
theory, which parameterizes the effects of small-scale corre- 
lations such as kinetic helicity on the evolution of the large- 
scale field (Moffatt 1978; Parker 1979; Krause & Riidler 
1980). Indeed, helicity has long been known to facilitate the 
production of large-scale fields under the condition that an ap- 
propriately defined dynamo number exceeds a critical value. 

Significant progress in the nonlinear formulation of large- 
scale dynamo theory has been possible by carefully design- 
ing and comparing simulations of turbulent magnetic fields 
with mean-field theory. It is important that both are applied 
to equivalent situations that are sufficiently simple (simple ge- 
ometries and boundary conditions, homogeneity of the turbu- 
lence, etc). Particularly illuminating is the case with triply- 
periodic boundary conditions where the energy of the large- 
scale field can exceed the energy of the small-scale field by a 
factor that is equal to the scale separation ratio, i.e. the scale 
of the system divided by the scale of the energy-carrying ed- 
dies. The phenomenology of this behavior based on magnetic 
helicit y conservation was already described by Brandenburg 
d2001l) . He found that full saturation occurs on a resistive time 
scale following a characteristic einschalt or switch-on curve 
pattern. The corresponding nonlinear mean-field the ory fo r 
helical dynamos was developed by Field & Blackman (2002), 
and extended to the case with shear by Blackman & Branden- 
burg (2002, hereafter BB02). 
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One of the important and surprising findings since the early 
works of Cattaneo & Vainshtein (1991), Vainshtein & Cat- 
taneo (1992), Gruzinov & Diamond (1994), and Bhattachar- 
jee & Yuan (1995) was the realization that the value of the 
magnetic Reynolds number, Re^, enters the nonlinear mean- 
field theory. Even today simulations have not yet been able 
to establish rigorously that the large-scale dynamo solutions 
obtained so far are asymptotically independent of the value of 
ReM- There are many aspects of this problem; the most impor- 
tant one is probably the possibility of catastrophic a quench- 
ing. This means that the field cannot saturate to equiparti- 
tion strength on a dynamical timescale and that in a closed 
or periodic domain fields of equipartition strength can only 
be reached on a resistive timescale. It is fairly clear now that 
catastrophic quenching can only be alleviated in the presence 
of magnetic helicity fluxes (Blackman & Field 2000; Klee- 
orin et al. 2000). What is less clear, however, is whether the 
helicity fluxes themselves depend on Re^ and on the mean 
field. 

Numerical evidence for large-scale fields in the presence 
of boundaries comes from simulations of both forced turbu- 
lence (Brandenburg 2005) as well as convection (Kapyla et 
al. 2008). The latter reference was particularly effective in ex- 
plaining the reason for the absence of a significant large-scale 
field in the simulations of Tobias et al. (2008), even though 
their boundary conditions would have allowed a helicity flux. 
The reason is that the helicity flux follows the direction of 
the contours of constant shear (Brandenburg & Subramanian 
2005b), but in the simulations of Tobias et al. (2008) these 
contours do not cross an open surface, because they used pe- 
riodic boundary conditions in the lateral direction. When us- 
ing instead open boundary conditions in the lateral direction, 
a strong large scale field is obtained (see Figure 17 in Kapyla 
et al. 2008). Alternatively, one can use vertical contours of 
constant shear, as was done in Kapyla et al. (2008). This was 
recently confirmed by Hughes & Proctor (2009) in an inde- 
pendent study. 

In the absence of magnetic helicity fluxes and without 
boundaries, strong large-scale fields can only be generated on 
a resistive timescale (Brandenburg 2001; Brandenburg et al. 
2001). With boundaries, and in the absence of magnetic helic- 
ity fluxes, the strength of large-scale magnetic fields decreases 
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like Re M with increasing values of Re^; see Brandenburg 
& Subramanian (2005b). Recent simulations of rigidly ro- 
tating convection have also been successful in generating a 
large-scale dynamo (Kapyla et al. 2009a). There, however, 
the ReM -coverage is not sufficient to determine a scaling for 
the saturation level. In the following we focus on the for- 
mer case using periodic boundary conditions, so no magnetic 
helicity can leave the domain, but large-scale fields can still 
emerge on a resistive timescale. We use such a configuration 
to study the quenching of the turbulent magnetic diffusivity, 
r/ t , and whether its value depends on Re^- The quenched 
value of T] t is crucial for determining the cycle period in an 
oscillatory <?Q dynamo. Indeed, BB02 showed that in the sat- 
urated regime and under the assumption of homogeneity the 
cycle frequency a> cyc of an aQ. dynamo can be written as 



-*cyc 



(i) 



where rji = ■% + rj is the total (turbulent plus micro- 
scopic) magnetic diffusivity in the quenched state and k m is 
the wavenumber of the mean magnetic field. Experiments 
with nonlinear mean-field models suggest that Equation (fl3 
remains approximately valid even in the mildly nonlinear 
regime. 

Equation ([TJ applies to homogeneous dynamo waves of suf- 
ficiently low amplitude, so that the field variation remains har- 
monic in space and time. Such behavior can only be mod- 
eled with periodic boundary conditions in the direction per- 
pendicular to the plane of the shear flow. The result is ex- 
pected to change in the presence of boundary conditions, but 
w C yc should still be proportional to r/j. Therefore, any Re^- 
dependent quenching of rjj would directly affect <w cyc , regard- 
less of whether or not magnetic helicity fluxes would alleviate 
catastrophic a quenching. This underlines the importance of 
studying the quenching of T] t . 

Determining the quenching of rj t has been attempted on var- 
ious occasions in the past. It is now generally accepted that, in 
two dimensions, rj t is catastrophically quenched proportional 
to 1 /ReM, provided the energy of the mean field is comparable 
to the energy of the turbulence. However, this is because in 
two dimensions the mean squared magnetic vector potential 
is conserved, which would not be applicable in three dimen- 
sions (Gruzinov & Diamond 1995). For three-dimensional 
turbulence, analytic theory predicts // t quenching proportional 
to the strength of the mean field independently of Re*/- Such 
a behavior also emerges from forced turbulence simulations 
where an initially sinusoidal magnetic field is found to decay 
proportional to exp(-?7T^m f ) (Yousef et al. 2003). The inferred 
quenching formula, applied to the corresponding mean-field 
dynamo models, is also found to produce the best fit to simu- 
lations that do show strong large-scale dynamo action (BB02). 
More recently, Brandenburg et al. (2008b) have computed the 
full aij and »/y tensors in the saturated state without shear, 
where the mean field is a Beltrami field. In that case, the tur- 
bulent magnetic diffusivity was found to be reduced by a fac- 
tor of about 5 as Re« is increased from 2 to 600. However, the 
case of Beltrami fields has certain limitations that are avoided 
when there is shear. 

2. THE SIMULATIONS 

We solve the stochastically forced isothermal hydromag- 
netic equations in a cubical domain of size (2jt) 3 in the pres- 
ence of a uniform shear flow, Uq - (0, S x, 0), with S = const, 

= -SA y x - (VU) T A + ifJ 2 A, (2) 



^ = -S UJ - c 2 s V lnp + -J x B + F visc + /, (3) 
T)t ' p 



£)lnp 

Dt 



-V -U, 



(4) 



where D/Dt — d/dt + {U + Uo) ■ V is the advective derivative 
with respect to the total flow velocity that also includes the 
shear flow, and [(Vi/) r A], = UjAjj in component form, F v ; sc 
is the viscous force, and / is the forcing term. As in earlier 
work (Brandenburg 2001), the forcing function is given by 



f(x, t) = Re{Nf k(t) exp[i*(f) • x + 10(f)]}, 



(5) 



where x is the position vector. The wavevector k(t) and the 
random phase -n < (p(t) < n change at every time step, so 
fix, t) is (5-correlated in time. For the time-integrated forcing 
function to be independent of the length of the time step 6t, 
the normalization factor N has to be proportional to 6r 1 ^ 2 . On 
dimensional grounds it is chosen to be N = foc s (\k\c s /6t) 1 ^ 2 , 
where fo is a nondimensional forcing amplitude. At each 
timestep we select randomly one of many possible wavevec- 
tors in a certain range around a given forcing wavenumber. 
The average wavenumber is referred to as kf. We force the 
system with transverse helical waves, 

fk = R • fi with R y = ; , (6) 
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where cr — 1 for the fully helical case with positive helicity of 
the forcing function, 



4 nohel) = (kxe)HV-{k-e)\ (7) 

is a non-helical forcing function, and e is an arbitrary unit 
vector not aligned with k; note that |/t| 2 = 1. 

We use periodic boundary conditions in the y- and z- 
directions and shearing-periodic boundary conditions in the 
^-direction. The main control parameters in our simulations 
are the magnetic Reynolds and Prandtl numbers, as well as 
the shear parameter, 



Re„ 



rjkf 



Pr 
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v -, Sh 
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By setting k\ = c s = po = fio = 1, we obtain dimensionless 
units of length, velocity, density, and magnetic field as 



[x] = fej , [u] = c s , [p] - po, [B] = V// poc s . (9) 

We solve the governing equations using the Pencil Code 3 
which is a high-order finite-difference code (sixth order in 
space and third order in time) for solving partial differential 
equations on massively parallel machines. 

The three essential parameters varied in the present study 
are the root-mean-square (rms) velocity (by changing the 
forcing amplitude fa), the relative helicity of the turbulence 
(by changing cr) and the microscopic magnetic diffusivity 
77. The rms velocity varies between 0.03 - 0.2, correspond- 
ing to /o = 0.01 - 0.05. The large-scale shear varies be- 
tween S = -0.05, . . . , -0.2 and leads to shear parameters 
Sh * -0.1, . . ., — 1.4. The high Sh cases studied in Branden- 
burg & Kapyla (2007) tend to produce very strong magnetic 
fields. The purpose of the study is to investigate the saturated 
state with different values of the mean magnetic field, B/B eq , 
so we need to regulate its value. One way to decrease the 

3 http: //www. nordita.org/software/pencil-code 
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Fig. 1. — J/ x , U y , B x , and B y at the periphery of the domain from a ran with Rcm s 209, Sh a -0.18, and Ptm = 10. Note the absence of any large-scale pattern 
in U y , but a rather pronounced pattern in B y with a fairly abrupt change of sign. 



magnetic fields is to increase the forcing amplitude to reduce 
the relative importance of the £2 effect. On the other hand, 
varying the value of <x can be used to regulate the a effect 
with given Sh and Re*?- This procedure was first adopted by 
Maron & Blackman (2002) and later by Brandenburg et al. 
(2002) in order to study the onset of large-scale dynamo ac- 
tion. In most cases <x is varied between 0.05 and 1, but small 
values do not always lead to clear large-scale dynamo action. 
Performing a large set of runs we have been able to cover a 
range of (B/B eq ) 2 from 0.3 to 270 and the magnetic Reynolds 
number varies between 1.5 and roughly 210. The numerical 
resolution of the simulations varies between 64 3 and 256 3 grid 
points. A summary of the runs is given in Table 1. 

3. FORMALISM 

In this section we briefly review some of the main predic- 
tions of the dynamical quenching model as applied to homo- 
geneous shear flows (BB02). In this approach, the quenching 
of a comes from an additional contribution that is propor- 
tional to the small-scale current helicity, which builds up at 
the same time as the a effect produces large-scale magnetic 
fields. The sign of the small-scale current helicity is opposite 
to that of the large-scale current helicity which, in turn, is also 
equal to the sign of the a effect. The evolution of the small- 
scale current helicity is calculated from the magnetic helicity 
equation. Using this equation together with the mean-field 



induction equation, BB02 derived expressions for the satura- 
tion amplitude, the cycle frequency, and the ratio of toroidal 
to poloidal field amplitudes in closed form. We begin with a 
more precise definition of k m in terms of the resulting mean 
magnetic field via 



k 2 m = fi {J ■ B)/(A ■ B), 



(10) 



where [Iq is the vacuum permeability, overbars denote mean 
quantities (later we shall specialize to horizontal averages), 
and angular brackets denote volume averaging. Preliminary 
results on turbulent diffusivity determined from the cycle fre- 
quency were presented in Brandenburg & Kapyla (2007). In 
the present study, we explore a much wider range of parame- 
ters and seek to understand the quenching behavior of rj t as a 
function of the magnetic field strength. 

For future reference we define at this point an analogously 
defined wavenumber of the fluctuating field via 



14=fio(J-b)/(a-b). 



(11) 



We also define effective wavenumbers of the fluctuating and 
mean fields via 



(12) 
(13) 

where €f and e m are the fractional helicities of the fluctuating 
and mean fields, respectively. We assume these two fractional 



e f k f = fioij ■ b)/(b 2 ), 
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helicities to be positive, but kf and k m can have either sign. In 
the stationary state the total current helicity must vanish for a 
closed system (Brandenburg 2001), which means that kf and 
k m must have opposite sign. Throughout this work the forcing 
function has positive helicity, so kf > 0, and therefore k m < 0. 

BB02 compared their re sults w ith those of direct simula- 
tions of Brandenburg et al. d2001l) . Theoretically, in the linear 
regime and at large magnetic Reynolds numbers, -k m can be 
as large as kf/2 (Brandenburg et al. 2002), but in the nonlin- 
ear regime, -k m will decrease until it reaches k\ (Brandenburg 
2001). Here we consider nonlinear solutions and will there- 
fore assume -k m = k\ . 

BB02 predicted that the saturation level of the mean mag- 
netic field is given by 



fin 
Bin 



ffkf 

£ m ki 



(1+ReT, 1 ). 



(14) 



Here, €fkf and e m k\ are the effective wavenumbers corre- 
sponding to the scale of the forcing and the mean field, re- 
spectively, and are defined via Equations dT2b and ( fT~3b . The 
equipartition value of the magnetic field is defined via 

Beq = W« 2 ) 1/2 - (15) 

According to the calculations of BB02, e m is directly propor- 
tional to the ratio of cycle frequency to shear rate. Further- 
more, e m is also proportional to the ratio of poloidal to toroidal 
magnetic field amplitudes. We have therefore multiple checks 
on the consistency of this simple model. 

4. RESULTS 

4.1. Saturation level of the magnetic field 

In most cases we find oscillatory solutions with dynamo 
waves propagating in the positive z direction. This is indeed 
expected from mean-field theory, according to which the di- 
rection of propagating is given by the sign of the product of 
a effect and shear. Indeed, positive helicity in the forcing 
should result in a negative a effect which, together with neg- 
ative shear, predicts a direction of propagation in the positive 
Z direction. 

Figure Q] shows the streamwise components of velocity and 
magnetic field from a run with Re^ m 209, Sh * -0.18, 
and Pi'm = 10. The velocity field is irregular while the y- 
component of the magnetic field exhibits clear large-scale 
structure. The x-component of the field also has a systematic 
large-scale component, but it is hard to see in a single snap- 
shot and without horizontal averaging because its amplitude 
is much lower than that of the y-component. 

This dynamo wave is well seen in animations showing the 
Z dependence of B y vs. time, but it becomes particularly clear 
when the field is averaged over the horizontal directions, indi- 
cated here by an overbar. An example of B x (z, t) and B y (z, t) is 
shown in Fig. [2] Note again the sharp sign changes of B Y (z, t). 
One can now see that the locations of these sharp sign changes 
coincide with the locations where B x (z, t) achieves positive or 
negative extrema. In fact, by comparing the two panels of 
Fig.|2]one can verify that a positive extremum of B x (z, t) leads 
to a change of sign of B v (z, t') from a positive value at ?' < t 
to a negative value at f > t, and vice versa. This is explained 
by the fact that shear is negative, i.e. S < 0, which turns a 
positive B x into a negative B y . 

The dynamo wave has a typical anharmonic shape (see 
Fig-EK just as has been seen before both in mean-field mod- 
els (Stix 1972) as well as in direct simulations (Brandenburg 



et al. 2001). Figure [3] shows flat positive or negative plateaus 
in By, during which the B x field was weak, but with a clear 

time derivative: while B y is positive, dB Jdt is also positive. 
A more careful look reveals that the plateaus are not com- 
pletely flat, but have a negative time derivative when B y is 
positive (and a positive time derivative when B x is negative). 
This relation can be interpreted as being due to a negative a 
effect in the relation 



dB 



V x (aB) + 



(16) 



Using the fact that there is a dynamo wave with propagation 
speed c, we have B = B(z - ct). Assuming constant a, this 
means that cB K = aB v . Since c is positive (compare Fig.|2]i, 

the opposite signs of the time derivatives of B x and B y sug- 
gest that a < 0. This is in agreement with the fact that the ki- 
netic helicity is positive and mean-field theory suggesting that 
for isotropic turbulence in the high-conductivity limit (Moffatt 
1978, Krause & Radler 1980) a is a negative multiple of the 
kinetic helicity. 

The oscillations are also discernible in the kinematic regime 
especially in runs where Re^ is small enough, so that the 
small-scale dynamo is not excited (see Fig . |4j) . The cycle fre- 
quency in the kinematic regime appears to be roughly con- 
stant in the range Re^ = 2 ... 10. In the following, however, 
we will be concerned with the nonlinear regime where the dy- 
namo has reached saturation. We study the dependence of the 
value of w cyc on Re M . Using Equation ([TJ we then calculate 
rjj and hence rj t . Figure [5] shows that the resulting value of rj t 
has a value roughly 2.5^0 for Re^ <: 4, where 



7to 



\tu 2 



(17) 



and r is the correlation time of the turbulence. The subscript 
zero refers to the kinematic case which is valid when the mag- 
netic field is weak. For small Re^ the correlation time is 
no longer determined by the turbulence but rather the micro- 
scopic diffusivity rj. Thus, in that limit r]i is expected to de- 
creases proportionally to rj. The markedly lower value of rj t 
for Re« ~ 2 can be interpreted in terms of this behavior. 

For Re^ » 1 the saturation formula given in Equation (TBI 
predicts that the ratio B^JB 2 ^ should be linearly proportional 
to €fkf/e m ki - (1 + ?7/?7to)- This is indeed in reasonable agree- 
ment with the numerical data; see Figure|6] 

BB02 gave two independent formulae for e m , one in terms 
of the resulting ratio of poloidal to toroidal mean fields, Q 1 , 

where Q 2 = <B?>/(^>, 



e m = V2Q" 1 = V2((l^>/<2^» 1 



n 



(18) 



and one in terms of the resulting ratio of the cycle frequency 
to the shear frequency, 

£ m = 2u> cyc /S. (19) 

Furthermore, assuming k m = -k\ in Equation (jT3j, we arrive 
at _ _ 

(J B) 



-Mo 



h{B 2 ) 



(20) 



We use the definition of e m given in Equation d20l > as the 
benchmark value to compare the other two expressions to. A 
scatter plot of the results is given in Figure [7] We find that 
the values given by Equations ( fT8l and ( fT9l are consistently 
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TABLE 1 

Summary of the different sets of runs. Here, 6 rms = Btms/Beq. (B x i = (B x )^ 2 /B^, and (B,,) 1 ' 2 = (B v )'' 2 /B, 



Set 


grid 


Re M 


Ptm 


klh 


Sh 




cr 




B 


ms 


ji2 
<B X ) 


1/2 




A 


64 3 


1.5 


1 


10 


-0.72 




0.2. . 


. 1 


8.0.. 


. 13.0 


0.14.. 


0.22 


7.2... 11.6 


B 


64 3 


46 ... 94 


10 


5 


-0.33 


0.16 


0.05. 


.1 


1.5. 


..2.4 


0.05 . . 


0.22 


0.95... 1.8 


C 


64 3 


27 ... 35 


10 


5 


-0.57...- 


0.44 


0.05 . . 


0.7 


2.0. 


..9.1 


0.05 . . 


0.09 


1.5... 8.2 


D 


64 3 


11. ..13 


10 


5 


-1.43... - 


1.16 


0.05 . 


.1 


2.7.. 


. 17.9 


0.03 . . 


0.10 


1.8... 16.2 


E 


128 3 


11 


10 


5 


-1.42... - 


1.36 


0.2. . 


. 1 


7.9.. 


. 14.4 


0.07 . . 


0.10 


7.2. ..13.1 


F 


128 3 


131... 209 


10 


5 


-0.29 


0.18 


0.05. 


.1 


1.2. 


..2.9 


0.04 . . 


0.20 


0.5... 2.2 


G 


256 3 


29 


25 


5 


-1.38...- 


1.29 


0.2. . 


. 1 


8.3.. 


.15.9 


0.05 . . 


0.06 


7.5... 14.3 


H 


64 3 . . . 128 3 


1.4. ..181 


0.5... 50 


5 


-0.13 ... - 


0.10 


1 




1.7. 


..3.1 


0.23 . . 


0.55 


1.2... 2.3 
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Fig. 2. — Space-time diagrams of B r (z, t) (upper panel) and B y (z, t) (lower 
panel). From a run with Re^ x 209, Sh » -0.18 and Ptm = 10. 
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Fig. 4. — Space-time diagrams of B x (z, t) (upper panel) and B y (z, t) (lower 
panel) from a run with Re/n ~ 1-4, Sh « -0.13 and Pr« = 1. 





B x (z=0,t)/B eq , 


B y (z= 


0,t)/B eq 


3 
2 
1 

if o 

-1 

-2 
-3 


* j \ / A V , u 

~~\"/,-^^'' \ \ 

- ---'si, \ J 






/ 'l * l / " 
/ A / - 






;''' Y " X S/ »- 




500 1000 


1500 


2000 


2500 3000 



Fig. 3. — Five times B x (dashed line) and B y (solid line) in the midplane, i.e. 
z = 0. From the same run as in Fig. [2] 



below the value of e m as given by Equation ( T20t . but that the 
trend as a function of cr is the same for all three quantities. 
Thus, we apply scaling factors of the order of unity for the 
expressions in Eqs. ( fT~8b and ( fl9] l in Figure|7]so that the data 
fall onto the same line. An exception are a few points for 
the cycle frequency, given by Equation (119t . where there are 
departures. The points deviating from the trend all belong to 
runs that have the weakest mean fields and occur for e m > 0.1. 

4.2. Effects from the vorticity dynamo 

A few runs exhibit signs of intermittent or continuous large- 
scale vorticity generation in the saturated state of the dynamo, 
reminiscent of a vorticity dynamo (Elperin et al. 2003; Yousef 
et al. 2008a,b). In Figure [8] we compare the evolution of the 
rms Mach number of the total (mean and fluctuating) veloc- 
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Fig. 5. — Cycle frequency (solid line) and turbulent diffusivity from Equa- 
tion {T} (dashed line) as functions of Re^ in the kinematic regime. 



ity both for hydrodynamic runs with and without shear, as 
well as a run with magnetic fields and shear. It is clear that 
in the absence of shear the rms Mach number reaches its fi- 
nal level quite quickly. During the first few tens of turnover 
times the same level is also maintained in the presence of 
shear, but the velocity gradually increases to much larger val- 
ues and reaches a new saturation level that is approximately 
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Fig. 6. — Bj^/B^ as a function of £fk{/e m k[ - (1 + fy/T/to)- The diagonal 
is shown for comparison. Note the reasonable agreement with theory (solid 
line). 
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Fig. 7. — Scatter plot of poloidal to toroidal field ratio (squares, scaled by 
factor 2.7) and normalized cycle frequency (plus signs, scaled by factor 2.4) 
vs. e ra . Note that for e m < 0. 1 both squares and plus signs scatter tightly 
around the diagonal. 



10 times larger. During that time a large-scale velocity field 
develops throughout the entire domain; see Fig.[9]for a typi- 
cal example of this behavior. In the hydromagnetic case the 
vorticity dynamo is quenched when the magnetic field grows 
to high enough level but in cases of weaker magnetic fields 
this quenching is often only partial as demonstrated by the 
dashed line in Figure [8] Unlike the anisotropic kinetic a ef- 
fect (Frisch et al. 1987), which is also suppressed by magnetic 
fields (Brandenburg & von Rekowski 2001), and works only 
for Re = Re^f/PiM < 8, the vorticity dynamo has been shown 
to work well up to Re m 100. A more thorough study of 
the vorticity dynamo can be found elsewhere (Kapyla et al. 
2009b). 

The overall behavior of such runs is quite different from 
those where a large-scale flow did not develop. We have there- 
fore excluded these runs from the plots for the sake of clarity. 
All in all the agreement between the three independent checks 
of e m is remarkably good especially for low values of e m . 
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Fig. 8. — Root mean square Mach number for two shearing runs with the 
same initial parameters (S = -0.2, v = 5 ■ 10~ 3 ) with (dashed line) and 
without (solid line) magnetic fields. The dotted line shows the result for a 
non-shearing hydrodynamic run for comparison. 
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Fig. 9. — U y at the periphery of the domain from a run with Rejf ~ 80, 
Sh ss -0.19, and Pr^ = 10, at a resolution of 64 3 mesh points at a time when 
the magnetic field is saturated and a large-scale vorticity dynamo develops. 



4.3. Turbulent diffusivity quenching 

According to the first order smoothing approximation 
(Moffatt 1978, Krause & Radler 1980), the turbulent diffu- 
sivity for isotropic turbulence in the high-conductivity limit 
is given by Equation ( fTTI ). Simulations of forced turbulence 
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without shear suggest that, independent of the value of the 
forcing wavenumber and amplitude, the Strouhal number is 
around unity (Brandenburg & Subramanian 2005a, 2007), i.e. 



St = k f u rms T as 1 



(21) 



In the kinematic regime, the magnetic difFusivity can therefore 
be written as 



J] to - 3 St u lms k f 



3 M rms^f 



(22) 



Recent numerical simulations employing the test-field proce- 
dure have shown that in the kinematic case, Equation ( 1221 is 
confirmed for Re^ between 1 and 200; see Sur et al. (2008). 
In order to study the quenching of i.e. the dependence on 
B, we normalize our results with the value of 7/ t o from a sim- 
ulation without magnetic fields. 

Our setup is similar to that us ed in a n umber of related sim- 
ulations by Brandenburg et al. (2008a) with non-helical tur- 
bulence, and Mitra et al. (2009) with helicity. In those sim- 
ulations the main focus was the determination of turbulent 
transport coefficients in the linear regime. In these studies 
it was found that the turbulent diffusivity increases by a factor 
of a few as |Sh| increases from 0.1 to unity. More recently, 
this work has been extended to the nonlinear regime, but so 
far only in the absence of shear (Brandenburg et al. 2008b). 
However, there is then a potential difficulty in that there could 
be additional terms in the functional form of the quenched a 
and j] t tensors that could mimic turbulent diffusion, so the split 
into two a and T] t coefficients is not unique. 

In the saturated state the growth rate of the large-scale field 
is zero. Thus, this situation corresponds to the marginally ex- 
cited state where the cycle frequency is given by Equation (1). 
Measuring therefore the dynamo frequency in the saturated 
state gives the (quenched) value of the turbulent diffusivity. 

Results for all of the runs are shown in Figure [10] The data 
scatter around the curve 



% = 



'7to 



(23) 



i+ g (\B\/B eq r 

with 7/ t o = (1.2-1.5) 77to, g = 0.3 and n = 1-2 being fit pa- 
rameters. A factor greater than one in the definition of fj t o re- 
flects the fact that in the kinematic regime, the turbulent mag- 
netic diffusivity is enhanced in the presence of shear (Bran- 
denburg et al. 2008a, Mitra et al. 2009). Note that the data 
for many of the runs with different magnetic Reynolds num- 
bers (ReM = 1 . . .210) seem to fall roughly on the same line. 
An asymptotic quenching of T] t inversely proportional to B in- 

—2 

stead of B has been predicted analytically by Kitchatinov et 
al. (1994) and Rogachevskii & Kleeorin (2001), and BB02 
found g ai 3 to be a good fit to the numerical simulations of 
Brandenburg et al. (2001). The fact that g is here different 
suggests that it is perhaps not a universal constant, but that it 
may depend on other parameters. 

The amount of scatter in Figure [10] does not depend in 
any systematic way on the value of Sh, except for the low- 
est values of Sh for which the values of j] t fall below the 
general trend. This, and the fact that rj t seems to reach val- 
ues above unity for low magnetic field strengths, might be a 
consequence of the turbulent diffusivity being enhanced when 
shear increases (Brandenburg et al. 2008a, Mitra et al. 2009). 
In the present case the unquenched value of M rms from a non- 
shearing run with the same forcing amplitude has been used in 
the definition of j] t Q. If the M lms from a hydrodynamic run with 
shear is used, the excitation of the vorticity dynamo would 
enhance its value by a factor of up to 5 (see Figure[8]l. 
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Fig. 10. — Turbulent diffusivity as a function of \B\/B e q. The squares and 
stars denote cases where -Sh < 0.2 (weak shear) and -Sh > 0.2 (strong 
shear). The upper curve is for n = 1 and i) t o/'7tO = 1-5, while the second 
curve is for n = 2 and i) t o/'7tO = 1-2, and g = 0.3 in both cases. 



In order to study the behavior of rj t as a function of Re^ a 
set of runs were performed where Pi~m was varied between 0.5 
and 50, keeping all the other parameters fixed, i.e. Sh « -0.1 
and Re = u lms /(vkf) as 4. The results are shown in Figure [TT] 
which demonstrates that rj t , as inferred from w cyc using Equa- 
tion ([1), decreases by a factor of about 2 when Re^ is in- 
creased from 10 to 180, while |B|/Z? eq varies from 1.7 to 2.5, 
with the larger values occurring at larger values of R&m- The 
low value of r\ x for the smallest value of Re^ might derive 
from the fact that the simple model used to estimate the value 
of the turbulent diffusivity may not be representative of this 
regime. On the other hand, there appears to be a declining 
trend of ?7t as a function of Re« approximately proportional 

1 /3 

to Re^ . However, higher Re^ simulations are needed to sub- 
stantiate this. On the other hand, data from nonlinear test-field 
calculations by Brandenburg et al. (2008b) show that 77, de- 
creases by a factor of 5 when Re M is increased from 2 to 600. 
Within error bars this is compatible with the present results. 

5. CONCLUSIONS 

Three dimensional direct numerical simulations of helically 
forced turbulence with imposed large-scale shear have been 
used to study the saturation level of the dynamo, its cycle fre- 
quency, and thereby the turbulent diffusivity and its quench- 
ing with the magnetic field. The parameters of the study were 
chosen such that a cyclic large-scale magnetic field develops 
with the dynamo wave propagating to the positive z-direction. 

We find that the saturation level of the energy of the mean 
magnetic field is compatible with a scaling of a quantity re- 
lated to the ratio of the fractional helicities of the small and 
large-scale fields, in accordance with Equation (TBI) and the 
prediction from nonlinear mean-field theory derived by BB02. 
Furthermore, the three independent measures of e m are in rea- 
sonably good agreement, especially for small values of e m . 

The turbulent diffusivity is found to be quenched approxi- 
mately inversely proportional to the large-scale magnetic field 
strength. This is in agreement with earlier analytical stud- 
ies (Kitchatinov et al. 1994; Rogachevskii & Kleeorin 2001). 
However, a small group of data points can also be fitted to a 
quadratic dependence on the inverse field strength. The de- 
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Fig. 11. — Turbulent magnetic diffusivity as a function of Re^, as obtained 
from the cycle frequency (solid line) for Sh = —0.1 and tr = 1, which trans- 
lates to ff » 1, compared with the corresponding result by Brandenburg et 
al. (2008b) using the test-field method (Sh = 0). The inset shows the Re« 
dependence of the cycle frequency <jj cyc (expected to be proportional to t/t; 
solid line) and the microscopic magnetic diffusivity (dotted line). 



pendence on Re^ is found to be weak which is also in ac- 
cordance with the analytical studies. More recent nonlinear 
test-field calculations (Brandenburg et al. 2008b), indicate a 
similar dependence. 

The present results may have implications for solar dynamo 
simulations. Clearly, the hope is that cyclic reversals occur on 
time scales that are asymptotically independent of the resis- 
tive time scale. What is less clear, however, is when asymp- 
totic behavior sets in. Evidently, even for the highest Re« 
simulations presented here the cycle frequency shows still a 
shallow decline with increasing magnetic Reynolds number; 
see the inset of Figure [TT] On the other hand, the correspond- 
ing change in the turbulent magnetic diffusivity is only by 



a factor of a few, even when the magnetic Reynolds num- 
ber changes by two orders of magnitude. This suggests that 
higher resolution simulations are needed to have any hope in 
seeing truly asymptotic behavior at large magnetic Reynolds 
numbers. 

A modest level of quenching of rj t would certainly still be 
compatible with solar dynamo models. In fact, already since 
the early 1970s it was clear that standard solar dynamo mod- 
els of aQ. type show cycle frequencies that are too high for 
realistic parameters, so the models predicted cycle periods of 
about byears instead of 22 years (Kohler 1970). A modest 
quenching of r] t would therefore point in the right direction, 
but one would hope that the decrease in j] t would eventually 
level off. Unfortunately, this is not evident from any of the 
simulations presented so far. 

As we have argued in the introduction, we do not expect 
the results for ^(Re^) to depend on the presence or absence 
of magnetic helicity fluxes. However, this expectation should 
be verified using simulations. Allowing for such fluxes would 
mean that we have to abandon the assumption of periodic 
boundary conditions. Although open boundary conditions 
would imply saturation on a dynamical time scale (not the 
resistive one, as in the present case of periodic boundaries) 
the saturation amplitude would scale inversely with Re^, if it 
was not for the shear-mediated helicity fluxes that allow for 
large saturation amplitudes on a dynamical time scale (Bran- 
denburg 2005, Kapyla et al. 2008). 
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